Tunnel splittings for one dimensional potential wells revisited 



Anupam Garg 

Department of Physics and Astronomy, Northwestern University, Evanston, Illinois 60208 

(February 1, 2008) 



Abstract 

The WKB and instanton answers for the tunnel splitting of the ground 
state in a symmetric double well potential are both reduced to an expression 
involving only the functionals of the potential, without the need for solving 
any auxilliary problems. This formula is applied to simple model problems. 
The prefactor for the splitting in the text book by Landau and Lifshitz is 
amended so as to apply to the ground and low lying excited states. 



1 



Typeset using REVTpjX 



I. INTRODUCTION 



The purpose of this article is to reexamine some formal aspects of the problem of calcu- 
lating the quantum mechanical tunnel splittings in a smooth, symmetric, one-dimensional 
double well potential, such as that in Fig. 1. Readers may justifiably wonder if anything 
new remains to be said on so mature a subject, and we assure them that, by and large, 
there isn't. The physical phenomenon is certainly very well understood, as are the basic 
mathematical ideas behind the calculations. We find, nevertheless, that especially in the 
case of the ground state splitting, there is some confusion regarding the correct WKB an- 
swer for this splitting, which is often taken in a form that is less accurate than it needs to 
be. In particular, we note that the formula for the splitting in the masterly text book by 
Landau and Lifshitz |IJ] leads to a prefactor which is incorrect for the ground state 0. A 
second problem is that the form in which this result is usually presented is poorly suited 
to calculation of the ground state splitting. In applying it to model problems such as the 
quartic double well potential, e.g., the unwary student is unnecessarlily led into asymptotic 
expansions of elliptic integrals, which must then be taken from standard tables of formulas. 
Because of this confusion, there is also confusion regarding the equivalence of the WKB 
and the instanton ||,|4| methods for calculating tunnel splittings. 

Like many other problems in physics of a similar nature, the above discrepancies and 
their resolution are at the same time "well-known" and not known at all! They are well- 
known to some experts. Thus, the correct answers for the splitting are implicit throughout 
early WKB studies of the anharmonic oscillator |5|||], and can be ferreted out with some 
work. They are also clearly known to many field theorists PHlOU, and to authors of more 



pedagogical articles pl]JT2|1 . At the same time, the confusion persists, and continues to crop 
up from time to time. Thus it seems worthwhile to address this issue here. 

It should also be stated at the outset that our article is of no interest if one only wants 
the splitting to "exponential accuracy", i.e., if one only wants the Gamow factor. In most 
physical applications, this is all that can sensibly be done, and the prefactor is best estimated 
as an "attempt frequency." From a mathematical point of view, however, the tunnel splitting 
is found as an asymptotic approximation in the limit h — * 0. An "exponentially accurate" 
answer for the splitting A is one for which one only has an asymptotically correct result for 
In A. An asymptotically correct answer for A itself requires worrying about the prefactor. 
More importantly, any formalism which did not give this prefactor, or which was incapable 
of giving it correctly even in principle, would be regarded as logically unsound. It is from this 
perspective that the prefactor is important. The confusion that is referred to above can be 
said to pertain solely to this prefactor, and readers who are not concerned with such arcana 
should stop reading here. On the other hand, the WKB analysis presented here is easily 
incorporated into a graduate level discussion of tunnel splittings, with little additional effort 
beyond that of the standard treatment, and there is no reason not to do so. The problem 
may also serve as a "real-life" physics example of some nontrivial asymptotic analysis. (To 
avoid misunderstanding, let us also state at the outset what we mean by WKB theory. We 
use the term in the broad sense used by Bender and Orszag [|HJ , or by Berry and Mount |14 



i.e., as a body of mathematical techniques that yield global understanding and systematic 
asymptotic approximations to solutions of many linear differential equations, 
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FIG. 1. Symmetric double well potential V{x) considered in this paper. The minima at x = ±a 
are taken to be quadratic. 

including Schrodinger's equation, and shares ideas with other methods such as asymptotic 
matching and patching, boundary layer theory, etc. Some readers may, however, define WKB 
to encomapss only the approximate exponential form ( |2.2| ) (and its oscillatory counterpart) 
and connection formulas at linear turning points. A subset of these readers may recognize 
that our treatment in Sec. II is tantamount to the use of quadratic turning point fT5| 



connection formulas (for which see Berry and Mount JT^] ) , and object that this goes beyond 
WKB. Such a semantic restriction of the scope of the term 'WKB', should in our view, be 
avoided. Since the method did not originate with Wentzel, Kramers, and Brillouin, but had 
antecedents in the work of Jeffreys, Rayleigh, Carlini, Green, Liouville, and perhaps others, 
there is no compelling historical reason for this restriction (as opposed to say, the usages 
'Einstein model' for lattice specific heat, or the 'Kronig- Penney model' for electrons in a 
periodic potential), and it is surely more useful to present the subject to students as one 
which allows for systematic improvement, and is not a closed subject of study even today.) 

It is with this motivation that we revisit the problem of calculating tunnel splittings in 
symmetric double well potentials. The article has three distinct aims. The first is to carefully 
calculate the prefactor multiplying the exponential of the turning-point-to-turning-point 
action integral in the standard WKB expression for the splitting A. The correct formula 
is Eq. ( |1.1| ) below. Our second aim is to present another formula [See Eqs. (0-0)] for 



the ground state splitting [jLq] , that reduces everything to the evaluation of two integrals 
involving the potential V(x), and does not require looking up any asymptotic expansions. 
The last aim is to show that the correct WKB and instanton MM methods do yield the 
same ground state splitting. We will do this by starting with Coleman's instanton method 
expression for the splitting and showing that it reduces to our simple formula. 

We expect that our WKB based disucssion will be accessible to students who have seen 
some graduate level quantum mechanics, even though our starting formula for the splitting 
is rarely ever mentioned in the common text books. One book where it does appear is 
again that of Landau and Lifshitz JlJ, who give an exceptionally lucid and self-contained 
discussion. We will comment further on this formula when we come to it in Sec. II, and give 
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a separate derivation of it in Appendix A. The instanton formalism, and the path integral 
approach to quantum mechanics upon which it is in turn based, are less likely to be familiar, 
but clear and accessible discussions have been given in this journal by Holstein [JTXJ] , and 



by Holstein and Swift |12| . The former discusses exactly the same problem as us, namely, 
the tunnel splitting, but for the special case of a quartic double well. We enthusiastically 
recommend Coleman's somewhat longer but authoratitive article on instantons |4[] to readers 
who wish to learn more about this technique. A related paper by Carlitz and Nicole JT7| 
may also be read profitably. 

A plan of our paper and summary of our results is as follows. We consider a potential 
which is smooth, reflection symmetric about x=0 [V(—x) = V(x)], and which has quadratic 
minima at x = ±a. The correct WKB answer for the ground state splitting is 



exp 



eir 



-a' n 



XV, 



Here, u is the frequency of small amplitude oscillations in the wells about x = ±a, ±a' are 
the classical turning points given by the equation 

V{a') = E = -hu + o{K), (1.2) 

and p is the momentum 

p(x) = [2m(V(x) -E)} 1/2 . (1.3) 

Note that p(x) is imaginary in the classically forbidden region —a' < x < a', and that we 
have taken V(±a) = 0. 

The proof of Eq. ( |1.1| ) is given in Sec. II. This entails matching the WKB wavefunction 
in the classically forbidden region to the exact harmonic oscillator wavefunction near the 
classical turning point, rather than use the connection formulas. This matching calculation 



may actually be found in a paper by Furry [19], but we include it because it is very short, 



and so as to have the complete argument in one place in a consistent notation. 

Prefactor corrections similar to those in Eq. fll.lj) accompany the excited states too, and 
are discussed in Appendix B. 

The expression ( |1.1| ) is not easy to use because the integrand in the exponential is close 
to a singularity near the limits, which means that the next to leading dependence of A on h 
is not manifest. In fact, the true preexponential factor varies as h 1 ^ 2 , and this is not obvious 
from Eq. (|1.1|) . In Sec. Ill, we shall, therefore, carefully extract the singular contributions 
to the action integral from the end points, and show that the splitting may be written for a 
general potential as 

/ 2\ 1 /2 

A = 2foJ^J e A e- s °/\ (1.4) 

where Sq (note the limits of integration) is the action integral 

S = f {2mV{x)) 1/2 dx, (1.5) 
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and 



.4 



muo 



2mV(x) 



x 



dx. (1.6) 



This expression does not have the complexities mentioned above. The h 1 ^ 2 dependence of 
the prefactor is apparent, and the answer involves only integral functionals of V(x). These 
formulas will be applied to two model problems in Sec. V. 

In Sec. IV, we turn to the instanton approach Here, the splitting is expressed as 

A = 2hK (H) 1/2 exp(-SoA), (1-7) 
where S is as in Eq. (|1.5|), and 



K 



det(-<9 T 2 + u 2 ) 1 1/2 
det'[-<92 + m-W'ixair))] 



;i.8) 



is the ratio of fluctuation determinants. Here, x c \(t) is the instanton, which obeys the 
(Euclidean) classical equation of motion 

~ V ' {Xd) = °' (L9) 

with the boundary conditions x c i(±oo) = ±a, and the additional condition x c i(0) = to fix 
the time translation degree of freedom. Further, V = dV/dx, V" = d 2 V/dx 2 , and the prime 
on the det in Eq. (|1.8|) means that the zero eigenvalue of the operator argument is to be 
excluded from a product of all eigenvalues. 

A short and pedagogical discussion of the instanton method may be found in Holstein 
|TT| . Holstein does not explain how to calculate the quantity K (which is proportional to his 
quantity Ki), so a few remarks that help elucidate its nature may not be out of place. Note 
first, that each determinant appearing in Eq. (L8) is of a one- dimensional Schrodinger-like 



operator, in which r plays the role of position, and either uj 2 or m~ 1 V"(x c \(T)) plays the 
role of the potential energy. Next we note, that the determinant of such an operator may 
be defined as an infinite product of all its eigenvalues. The spectrum of operators at hand 
may be rendered completely discrete by adding hard walls at r = ±T, and letting T — > oo 
at the end. This scheme has the advantage that it allows us to put the eigenvalues of both 
operators in one-to-one correspondence. This in turn enables us to argue that one eigenvalue 
is missing in the denominator, and since the eigenvalues have the dimensions of d 2 or uj 2 , 
we see that K has the dimensions of frequency. For a smooth potential V(x), the curvature 
m^ 1 V ; '(x c [(r)) is never very different from uj 2 in the interval —a < x < a, so there is only 
one frequency scale in the problem. It follows that K is of order uj. 

To the author's knowledge, the equivalence of Eqs. (|i~T| ) and (|1.7|) has only been shown 
for particular examples, such as the quartic double well potential. A general demonstration 
is lacking, although it is implicit in the fact that both methods start from the same point 
and are correctly executed. A direct demonstration is therefore of some value. Secondly, 
Eq. ( |1.7| ) is complicated just as Eq. ( |1.1|) was. Indeed, the ratio of determinants seems 
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more intimidating and harder to calculate than the integral in Eq. ( |1 . 1|) . Langer's original 
calculation || for K is quite involved, and while it has been greatly simplified by Coleman 

it still requires solving the auxilliary problem of finding the instanton x c \{t). 

To relate the instanton and WKB answers, ws therefore show in Sec. IV, that Coleman's 
result for K may be expressed in terms of the integral A in Eq. ( |1.6| ). Thus, although we 
do not provide a complete derivation of the instanton result for A, once Coleman's formulas 
are accepted, our analysis shows the equivalence of the two approaches. 



II. WKB FORMULA FOR GROUND STATE SPLITTING 



Our starting point is the formula [20 



2h 2 

A = —^,(0)^(0), (2.1) 

777 

where if)o(x) is an approximate solution to Schrodinger's equation with energy Eq, which 
is localized in the right hand well, and decays away from that well, in the entire central, 
classically forbidden region. Further, it is normalized to yield unit total probability in the 
right hand well. Note that it is not necessary to examine or specify the behaviour of i[>o(x) 
near x = —a'. Lastly, ip' (x) = dipo(x)/dx. 

We digress briefly at this point to comment on the formula ( |2.1| ) , which does not seem to 
be widely known. It is sometimes named after Conyers Herring, who derived it (in a slightly 
more general form, in fact) in the course of evaluating the g-u splitting of the two lowest 



electronic states of the H2 molecular ion [21,22]. The actual derivation is simple, and very 



similar to the usual half-page of analysis used to argue that the the Hamiltonian operator 
in the position representation is Hermitean and has real eigenvalues, and so we refer readers 
to Landau and Lifshitz [l[] for the details. Instead, we present an alternative derivation in 
Appendix A. 

Resuming our argument, we note that the WKB approximation for ipo{x) in the region 
(a' — x) > (Ti/muj) 1 / 2 is 

C 1 f x 

4>o{x) = ^=J= exp - \p(x')\dx', (2.2) 



v(x)\ 

where v(x) = p(x)/m, and Co is a constant to be found by matching on to the solution in 
the well. To the accuracy of this solution, V , o( x ) ~ (m\v(x) \/h)ipo(x), so that 

A = 2hC 2 . (2.3) 

To find Co, we first note that near x = a, ipo(x) is very accurately given by the ground 
state harmonic oscillator wave function 



^o(x)=( v — ) exp 



muj 2 
•— -{x — a) 



(2.4) 



In fact, this form holds well into the forbidden region, and so can be directly matched to 
Eq. (|2.2|) without invoking connection formulas. In the overlap region, we may write 
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\p(x)\ = muj (a — x) 2 — 



1/2 



(2.5) 



where 



(2.6) 



The term Mq in Eq. ( |2.5| ) may be neglected in evaluating the \v(x)\ 1//2 prefactor in Eq. fl2.2|) , 
so that we may write 



Co 



uj(a — x) 



exp 



1 

h Jo 



\p{x')\dx' + <&{x) 



(2.7) 



where 



(a — x ) 



i\2 



Ur 



1/2 



dx 



muj(a — x) 2 , f 2(a — x) 
- + In -i ^ 



1/2 



Comparing Eqs. 



2^ \ Mo 

and (O) with Eq. ([O]), we obtain 



4 \a — x 



(2.8) 



Co 



47re 



C 4 


■ I ra' 


exp 





(2.9) 



Substituting Eq. ( p.9|) in Eq. ( |2.3|) , and making use of the reflection symmetry of V(x), 
we obtain 



A = — , exp 



a ' \p[ 

-a' k 



dx 



(2.10) 



which is Eq. 

Landau and Lifshitz obtain Co by using connection formulas near a linear turning point, 
and the standard WKB result for the normalization of a bound state. For the ground state, 
this procedure is not accurate enough, and we must proceed as above. (Some readers may 
recognize in our procedure, the use of quadratic turning point formulas Once this is 

realized, similar prefactor corrections are expected for the low lying excited states. These 
corrections are obtained in Appendix B. 



III. SIMPLER EXPRESSION FOR GROUND STATE SPLITTING 

In this section we will show that the WKB expression leads to the simpler result 



( |1.4| - |1.6| ). The objective, clearly, is to let the action integral run from —a to a instead of 
from —a' to a'. To that end, let us denote y = a — x, V(x) = U(y), and define, for a general 
u, 



/(«)= [ [2m(U(y)-U(u))] 1/2 dy. 



(3.1) 
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The quantity I(uq) is clearly half the action integral in Eq. (|1.1|) , i.e., 

I(u ) = ^£ ai \p(x)\dx. (3.2) 

Our goal is to expand I{u) for small u. The analysis in Sec. II reveals [see Eq. Q2.8| )] that 
this expansion contains a term varying as u 2 lnu, which means that the expansion cannot 
be found by simply differentiating I(u). More precisely, the term of order u 2 lnu is easily 
found, but the term of order u 2 , which we also need, is rather harder to get. 
We therefore resort to a subtraction and write I(u) = I\(u) + I[(u), where 

h(u)= [ a yj2mU(y)dy, (3.3) 

I'M = £(^/2m[U(y) - U{u)\ - ^2m£%))dy. (3.4) 

I\{u) can be directly expanded: 

h{u) « £ ^2mU(y)dy - £ (muoy + 0{y 2 )) dy 

= 1(0) - ^muu 2 + 0(u 3 ), (3.5) 

while for I[(u) we have, as u —> 0, 



a 



I'Au) w -{muouf I " jL= . (3.6) 

^ ^2m[U(y) - U(u)\ + s]2mU(y) 

Note that in writing this expression, we have taken U(u) ~ muj 2 u 2 /2 in the numerator, as u 
is small. It still can not be expanded directly, however, so we perform another subtraction, 
and write I[ — I2 + I3, where 

IM = -Umujuf r dy . (3.7) 



2 



y/2m[U(y) - U(u)} 



The difference J3 contains a proportionality factor U(u), which we again approximate by 
rruj 2 u 2 /2. This leads to 



h(u) ~ — (muJu) A I J 



2 



u pm[U{y) - U(u)] y2m[U(y) - U(u)] + pmU{y) 

Let us consider ^(u) first. We can clearly write the integral as 
dy 



2 ■ 



(3.8) 



u mujy/y 2 — u 2 



+ 



^2m[U(y) - U(u)} mu^y' 2 - u 1 



dy. (3.9) 



The first integral can be done exactly, while in the second we can set u = inside the 
integrand and in the limits to leading order. We thus find 
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Ioiu) ~ — mijju 2 hi (muu) 

2K ' 2 u 2 



^2mU(y) muJ V 



dy + 0(u 4 lnu). (3.10) 



The leading behaviour of 13(11), on the other hand, is controlled by the lower limit in the 
integral ( |3.8| ), where we can again approximate U(y) by muj 2 y 2 /2. Hence, 



h{u) ~ -muu 



dy 



u ^Jy' 1 — u l \Jy 2 — u' 2 + y 



-muju 2 + 0(w 4 ). 
4 



(3.11) 



Collecting together Eqs. Q3.5Q , ( |3.10|) , and ( |3.11|) , and putting u = u = (/j/mw) 1 / 2 , we 



get 



, , . . h , 2a h h f a 
J(uo)«/(0)--1b 



2 m 4 2 Jo U2mU(y) V 



muo 



dy + 



(3.12) 



The last integral in the above equation is nothing but the quantity A defined in Eq. ( |L6|) . 
Using m = (h/muj) 1 / 2 once again, and putting together Eqs. ( |3.12|) , fl3.2|) , and (|1.1|) , we get 



A 



2huo ( muia 



2\ 1/2 



Since 1(0) = S /2, this is nothing but Eq. Q). 



\/ei 



A_-2/(0)/» 



ee e 



(3.13) 



IV. EQUIVALENCE OF WKB AND INSTANTON RESULTS 

Our next step is to show that the instanton result (|1.7|) also leads to Eqs. (|1.4H1-6|) , and 
thus prove the equivalence of the WKB and instanton results for the ground state splitting. 
To do this, we use Coleman's result for the ratio K H. According to him, 



K = V2u(3, (4.1) 
where (3 is related to the asymptotic, r — > ±00 behaviour of the instanton velocity via 

Xl (r) = (^) 1/2 ^^(3e^ T asr^ioo. (4.2) 
V So J dr 

It is easy to integrate the equation of motion ( |1.9| ) for the instanton, and obtain X\(r). 
Using the fact that x c i(0) = 0, we have 

pci dx 

r = m (4.3) 
l2mV(x) 



This diverges as x c \ — > a, as it should. We extract the divergence by subtracting and adding 
the leading singular part of the integrand. This yields, 
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t = m 



^2mV(x) muj(a-x) 



, 1, / a 

ax H — m 



lo \ a — x c \ , 



A 
+ - 



I f a 
— In 



as x c \ a. 



A is, of course, the quantity defined in Eq. (|I~6|). 
Thus, as r — ► oo, 



£ci( r ) ~ a — ae^e WT , 



cl A — u;r 

« owe e . 



Comparing with Eq. (|4.2|) , we can read off /3 immediately: 

(3 = au ^— ) e A . 

Hence, 



and 



if = acj ( — - * ! 

D 



A = 2hu (^^e^-^, 



which is what we set out to show. 



V. ILLUSTRATIVE EXAMPLES 

We conclude with two elementary examples to which we apply Eqs. ( |1.4| - |1.6| ) 
The first example is that of the quartic double well, 

V(x) = V (x 2 - a 2 ) 2 /a\ 

Here, Vq is the barrier height. The frequency uo is given by 

co = (8V /ma 2 ) 1/2 . 

It is simple to perform the integrals, and show that 

S 4 / '2mV a 2 \ 1/2 16 V 



h 3 v h 2 

A = ln2. 



3 hoo' 



Further, {mua 2 /nh) 1 / 2 = (3S /27ih) 1 / 2 , so that 
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/ C x 1/2 

\2iih) 

which is a well known form for the answer. 

Our second example involves the spin Hamiltonian 



(5.5) 



H = -jS 2 z - aS x , (5.6) 

where S a (a = x, y, z) are components of the dimensionless spin operator S obeying the 
commutation rules [S a , Sp] = ie a pySy, and we also take a > 0, 7 > 0. 

In the limit where the magnitude S of the spin is very large [23] , Eq. ( |5.6| ) can be viewed 
as tending to a classical Hamiltonian, wherein the dynamics are defined by giving the Poisson 
brackets [S a , S^pr = Ca/Zh/Sj- The classical energy has two degenerate minima when the spin 
lies in the xz plane at angles 9 oi n — 9 to the z axis, where 

sm6 = a/2~fS. (5.7) 

For the quantum mechanical problem, with large but finite S, we expect the ground states 
around these classical orientations to be admixed by tunneling, thus giving rise to a small 
splitting of the energy levels. 

Compared to the tunneling of massive particles with a position coordinate, the tunneling 
of spins is somewhat less familiar, but it is a perfectly bona fide instance of the general 
tunneling phenomenon. Spin tunnel splittings have been calculated by a variety of means for 
some time now p4[^5j, but obtaining the prefactors in the tunnel splittings correctly to order 



S° as S — > 00, has proven to be a fairly difficult task, and the instanton calculations fl24|.|2l 
are especially subtle. The discrete WKB method uses only elementary methods of analysis, 



but the calculations to date [0,0 are lengthy. There is an even simpler method, 



however, which is due to Scharf, Wreszinski, and van Hemmen ||28|| . These authors expand 
a general state in the 5^ eigenbasis (S z \m) = m\m)) as 1^) = J2 m D m \m), and 

construct a generating function J2 m D m x m . After a few changes of variables, the Schrodinger 
equation for is turned into the following Schrodinger equation [See their Eq. (2.17)] for 
a wavefunction y(z) related to the generating function: 



l 2 ^ + V{z)y{z)= 1 Ey{z) 1 (5.8) 



with the potential 



V(z) = -ct 2 (cosh z — cosh zq) 2 , (5.9) 

coshzo = (25+ 1)7/0. (5.10) 

The quantity E in Eq. ( |5.8p is the energy eigenvalue of the Hamiltonian ( p.6| ). 

The potential (|5.9| ) is even about z — 0, and has minima at ±Zq. Our formalism is 
directly applicable if we identify h 2 /2m with 7 2 , and muj 2 with V"(±z Q ) = ±a 2 smh 2 z . The 
action integral Q1.5Q is easily seen to be 

Sn a r z o 



Dq Oi f z ° 

— — — I (coshzo _ cosh z)dz 
ft 7 JO 



7 

OL 

— (z coshz — sinhzo), (5.11) 

7 
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while the correction (|1.6|) is given by 



.4 



-0 



sinh Zq 



cosh zq — cosh z 



Zq- Z 



dz. 



(5.12) 



Since the integrand is nonsingular at z = Zq (or if one so wishes to view matters, has a 
removable singularity at that point), the integral can be found by changing the upper limit 
to zq — 5, integrating the two terms separately, and letting 5 — > at the end. The integrals 
involved are elementary, and the final result is 



A = In 



2 sinh Zq * 



(5.13) 



It remains to substitute Eqs. (|5.11| ) and (|5.13|) into our general formula (|1.4j ). Recalling 



the factor of 7 on the right hand side of Eq. ( |5.8| ), and the equivalences h /2m 



muo 



V"(zq), we obtain 



4a 3 / 2 sinh 5 / 2 



V2t^ 



exp 



a 



(zq cosh Zq — sinh z ) 

7 



r 



(5.14) 



It is easier to interpret this result if it is cast in terms of the angle #0 in Eq. (577). After a 
certain amount of tedious algebra, one obtains, correct to order S° as S — ► 00, 



A 



87s 3 / 2 cos 5 / 2 6 / 1-cosflc/ 
ttV2 S i n Qq [ 1 + cos# 



5+1/2 



,2S cos f 



(5.15) 



This is precisely the form quoted in Ref . p7j • Its relation to previous calculations is discussed 
there. 

From the viewpoint of this article, the interesting comparison is with Scharf, Wreszinski, 
and van Hemmen [B?|. These authors only consider the case where a and 7 are fixed as 
00, so that 6*o = 0{1/S). In that limit, they follow the Landau-Lifshitz prescription, 



S 



and obtain a ground state splitting A' = (2S/n)e J , with 
J = -2S In (--f 2 S 2 ) +2S + {2S + 1) In f 7 (25 + \ 



- - In 

2 



7 2 (S + i) 



(5.16) 



With a little work one can show that A'/A is exactly (e/ir) 1 ^ 2 . The discrepancy can in fact 
be seen in the comparison between numerical and analytical results in Table 1 of Ref. PS| . 
While (e/ff) 1 / 2 = 0.930, the ratio of their analytical answer for the splitting to the numerical 
one (see the last two columns) decreases from 0.987 to 0.951 as S increases from 5 to 11, 



which is quite close. A Richardson transformation [[29|] of this ratio does suggest that it 
tending to 0.93, but one cannot be certain of this conclusion, as there are two erratic values 
at 0.90 and 0.91. 
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APPENDIX A: HERRING'S FORMULA 



Let us denote the two states localized in the left and right wells by \R) and \L) respec- 
tively. These states are degenerate in the absence of tunneling with an energy Eq. The 
tunnel splitting is A. The Hamiltonian matrix in this two-state subspace is given by 

«-{!hf) < A1 > 

Assuming that the system starts in the state \R) at t — 0, it is straightforward to show that 
the probability -Pr(£) for finding it in the same state at a later time t is given by 

P R (t) = cos 2 (At/2h). (A2) 

In particular, 

f^A^y (A3) 

To relate this abstract space description to that in position space, we make use of the con- 
tinuity equation for probability. For a general wave function i[>(x, t) that obeys Schrodinger's 
equation, the probability density P(x,t) = \if)(x, t)\ 2 obeys 

T = -^ Im ^- f ™< < A4 > 

where ip' = dip/dx. The probability -Pr(£) for being in the right well is then given by 

POO 

P R (t) = / P(x,t)dx. (A5) 
Jo 

Differentiating with respect to t, making use of the continuity equation, and the fact that 
ip(x, t) — > as x — > oo for any well behaved solution, we obtain 

^ = ^lm[f(0,^'(0,f)]. (A6) 

Let us finally consider the states \R) and \L) in position space. It is evident that we should 
take (x\R) to be the function described as ipo( x ) i n Sec. II. By symmetry, (x\L) = iPq(-x), 
and the energy eigenstates are (ipo(x) ± x))/V2. It then follows that with the same 
initial conditions as above, the wave function at an arbitrary time t is given by 

if)(x, t) = ip (x) cos(At/2h) + iip (-x) sin(At/2^). (A7) 

Hence, 

Im[<0*(O, t)ip'(0, t)} = ^ (0)^o(0) sin(At/^). (A8) 
Substituting this result in Eq. ( |A6| ) and comparing with Eq. (^), we obtain 

A = (2ft 2 /m)Vo(0)Vo(0), (A9) 
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which is Eq. (|2.1| ), Herring's formula. 

It is apparent from our derivation that Herring's formula holds whenever the energy 
eigenf unctions are well approximated by the combinations (ipo( x ) =t V ; o( — a; ))/v / 2- (More 
precise statements of the conditions for its validity are given in Ref. IffiH .) Use of this formula 
greatly simplifies the labour required to solve all the standard double-well problems: the 



double square well with infinite side walls [pO[pT| , the double delta function |31| , and the 



double harmonic oscillator with a kink [3Hfl. For smooth potentials, where WKB is indicated, 
it is far superior to the approach where one uses connection formulas at all four turning points 

. Further, the formula brings out the physical point that a tunnel splitting is intimately 
related to a tunneling amplitude, since it relates the amplitude to make a transition between 
wells per unit time, —iA/h, to the probability current. 

The reader will also have noticed that the above argument cannot be carried out for an 
asymmetric potential, and indeed the very concept of tunnel splitting is then very delicate. 
A tunneling amplitude can of course still be defined, but since this amplitude is exponen- 
tially small in general (on account of the Gamow factor), mixing between left and right 
well states will be negligible unless the bottoms of the two wells are tuned to the same 
exponential sensitivity. The mathematical formulation of these points leads to extremely 
unpleasant transcendental equations, and the situation is not significantly improved in the 
case of symmetric potentials if one approaches the problem solely in terms of enforcing con- 
tinuity of the wave function and its derivative. It is probably because of this fact that most 
introductory or intermediate quantum mechanics texts do not consider double-well tunnel 
splittings when they discuss the WKB method. Use of Herring's formula (with or without 
WKB) would make the problem much more tractable, and its widespread adoption is thus 
greatly to be desired. 



APPENDIX B: WKB FORMULA FOR HIGHER STATE SPLITTINGS 

It is to be expected that the formula given by Landau and Lifshitz PJ is increasingly 
accurate for higher pairs of states, assuming of course, that the WKB approximation is still 
applicable. We will find that this is indeed so. The splitting for the nth pair of levels, A n , 
is given by 



Two 

A n = g n — exp 

7T 



-^-dx 
n 



(Bl) 



where ±a' n are the classical turning points for the nth energy level pair (n + h)hu, and 



9n 



nl 



( n+ iy 



"2 e -(n+i)_ 



(B2) 



Note that if g n were unity, we would have the formula of Ref. [0 . The corrections are indeed 
small: go = (ix/e) 1 ! 2 ~ 1.075, g\ ~ 1.028, g^ ~ 1.017, and so on. Stirling's formula for n\ 
shows that g n — > 1 as n — > oo. 

The derivation of Eq. ( |B1| ) and the conditions under which it holds are straightforward 
though long. The starting point is still Herring's formula, Eq. (|2.1| ). The procedure of Sec. II 
can be applied word for word, with a', uq, and Co replaced by a' n , u n , and C n , respectively. 
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The essential part is to find ip(x) in the classically forbidden region and match it onto the 
nth harmonic oscillator state in the vicinity of the well x ~ a. We leave it as an exercise 
- or see Ref. [JHJ - to show that the function defined in Eq. ( |2.7| ), is given in the 

overlap region by 

. Imuj . , 2 2n+l f 2(a - x)\ 2n + 1 
•(«) = - 2 — (« - «) 2 + — 1" ("4^ j + — + - ■ < B3 > 

so that the leading approximation to i/jq(x) is given by C' n (a — x) n exp(— x 2 /2uo), where 

C n / 2 \ n+ i (I r< 2n + l\ 

This is of the precise form required to be match on to the nth excited state harmonic 
oscillator wave function: 

with £ = (a — x)/uq. The match yields the constant C n , and the splitting A n , which equals 
2h\C n \ 2 , is easily shown to be given by Eq. ( |B1| ) and (|B~2j). 

The calculation sketched above assumes that the turning point and the wavefunction 
in its vicinity are well approximated by taking the potential well to be parabolic. This is 
clearly less accurate as n gets large. Defining y and U (y) as in Sec. IV, let us keep the cubic 
and quartic terms in y in U (y) : 

U(y) = ^muj 2 y 2 + ay 3 + {3y\ (B6) 

The wavefunction for the nth state can be found by perturbation theory assuming that a 
and (3 are small. The key requirement is that the correction terms be small compared to 
the unperturbed wavefunction near the turning point y = u n . The dominant corrections are 
those that entail the Hermite polynomials of order n + 1 through n + 4. Since the coefficient 
of the largest power of x in these polynomials is known, the conditions for the corrections 
to be small are not difficult to find. They are 

e*WTmrrw¥-- < B8) 

(6n + 7)(2n + 1) n 

It is also possible to show that when these conditions hold, the corrections to u n and the 
energies of the harmonic oscillator states are negligible, so that the conditions are self- 
consistently derived. 



It need hardly be said that the conditions (p7| ) and ( pq) are only necessary, not sufficient, 



for Eq. ( [B 1| ) to apply. In addition, one must have the requirement for the WKB approxima- 
tion itself to hold, which in the present context can be stated as (n + \)hu V(0) — V(a). 
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